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' Abstract In this paper, the spectra of the persistent emission from anoma- 

' lous X-ray pulsars (AXPs) and their variation with spin-down rate is 

, considered. Firstly, based on an accretion-powered model, the influences of 

• both magnetic field and mass accretion rate on the spectra properties of 

^ ' AXPs are addressed. Subsequently, the relation between the spectral prop- 

|H I erty of AXPs and mass accretion rate M is investigated. The result shows 

, that there exists a linear correlation between the photon index and mass 

^ ■ accretion rate, and the spectral hardness increases with increasing AI. A 

' possible emission mechanism for the explanation of spectral properties of 
H ' 

' AXPs is also discussed. 
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1 INTRODUCTION 

Anomalous X-ray pulsars (AXPs) are one of the enigmatic class of Galactic high energy 
sources. These sources differ from known magnetic accreting X-ray pulsars found in high 
and low mass X-ray binaries (HMXBs and LMXBs) (Mereghetti & Stella 1995; Van 
Paradijs et al. 1995). AXPs are sources of pulsed X-ray emission with spin periods in the 
6-12s range, very soft X-ray spectra, secular spin down on time scales of ^ 10"^ — 10^ yr, 
and lack of bright optical counterparts. Two or possibly three of AXPs are close associated 
with supernova remnants (SNRs). Additionally, AXPs share some similarities with the 

Soft Gamma-ray Repeaters (SGRs) (Hurley 2000). 
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Two broad classes of models have been proposed to understand the X-ray emission 
from AXPs. The first class of models assume that the sources of AXPs are isolated 
neutron stars with ultra-magnetic field strengths in the range 10^'* — lO^^Gauss - i.e. 
"magnetars" (Thompson & Duncan 1996; Hely & Hernquist 1997). The spin-down of 
the pulsars is primarily due to magnetic dipole radiation, then the AXPs have enor- 
mous surface magnetic dipolar fields. The second class of models suggested to explain 
AXP emissions is that they are accreting from a disk of material leftover from the su- 
pernova explosion that created the neutron star (Chatterjee ct al. 2000, hereafter CHN; 
Alpar2001; Lu & Cheng 2002; and Marsdcn et al. 2001). Such models do not require 
neutron stars with unusually strong magnetic fields. In these cases, the spin-down torque 
is external and, supposedly, a neutral consequence of accretion braking (Francischell & 
Wijers 2003). A very different scenario, based on strange matter stars, has been proposed 
by Dar & DeRujula (2000). They argued that AXPs are either strange stars or quark 
stars in which the X-rays are powered by gravitational contraction, and the spin-down 
is due to the emission of relativistic jets. 

AXPs could be magnetars, but valid alternative models cannot be ruled out with a 
decisive proof (Mereghetti et al. 2002). In the magnetar model, a more stable luminosity, 
spectrum and P noise can be expected more naturally than in accretion-powered pulsars. 
The spectra and light curves expected from the surface of highly magnetized neutron stars 
have been computed by several authors (Ozel 2001; Zane et al. 2001; Ho & Lai 2001). The 
spin-down of magnetars is due to a combination of standard magnetic dipole radiation 
torque and torque from the wind (Harding et al. 1999). Both of these torques increase 
strongly with magnetic field strength, and therefore the hardening of the power-law 
spectral component with O implies a similar hardening of the underlying Alfvcn wave 
spectrum with increasing magnetic field B. It is known that the observation of AXPs 
shows that the spectral hardness decreases with increasing spin- down rate (Marsden & 
White 2001). Therefore, the correlation of the spectral hardness (photon index T) and 
the spin-down rate (Cl) predicted by magnetar model is inconsistent with the data. 

Based on the above investigations, alternatively, we relate accreting neutron stars 
depending on both the magnetic field and the accretion rate M to model the spectral 
properties and spin-down rates of AXPs. The purpose of the present paper is twofold. 
Firstly, based on the observation, the relation between the photon index F of AXPs/SGRs 
and the accretion rate M is investigated. Secondly, a possible emission mechanism is 
suggested to model the spectral properties of AXPs/SGRs. Finally, some predictions 
based on the relation between T — M are discussed. 



anomalous X-ray pulsars 3 

2 PARAMETERS AND OBSERVATIONAL DATA FOR AXPS 

An important parameters in describing AXPs is its spin period (P). It is known that 
the AXPs slowly brake, so the time derivative of the period, P, is another one. In ad- 
dition, based on two observations, other useful pulsar parameters can be defined: the 
"characteristic age" t = P/2P, and the magnetic field 

where / ~ lO^^gcm^, and r ^ 10® cm. 

The role of period derivatives P is confirmed when observing the AXPs spectral index 
variation. Following the data analysis of Marsden & White ( 2001) for the phase-averaged 
spectra of the AXPs and SGRs in the range of 0.5 — 10.0 keV, we obtained the basic 
data of all SGRs and AXPs with known spin periods and period derivatives detected 
with ASCA which is summarized in Tabled The period derivative of one object AXP 
1048-59 cannot be determined definitely, so we take two values throughout the following 
analysis. And the photon indexes here are derived in terms of two-component black-body 
plus power law spectral models. Figure^ shows the variation of the photon index F vs. 
spin-down rate VL for such AXPs and SGRs. 

The relation of the photon index F and P, P is fitted in terms of the linear equation 
F = alogP -f /31ogP + 7. For five AXPs and two SGRs listed in TableQ the fitting of 
F(P, P) is 

F(P, P) = 2 log P - log P12 + 2.5 , (2) 

where P12 denotes P/lO'^^ss-^ Substituting Q. = 27r/P and tl = 2tiP/P'^ into Eq. El 
the correlation between F and Cl can be obtained. F is a linear function of Vl. The photon 
index F decreases with increasing spin-down-rate f2. (see Fig. 

Furthermore, the relation ((T — B)) between the photon index and magnetic field 
strengths for such AXPs/SGRs is studied. We plot it in Figure Figure 2 shows that 
the harder the spectra of AXPs/SGRs, the stronger the magnetic fields. It seems that 
the present magnetar model cannot give this result yet, so it may require the further 
theoretical considerations. The following study shows an accretion-based model (CHN) 
could give a possible explanation to the observed correlation between the photon index 
and spin-down rate in AXPs and SGRs. 

3 THE SPECTRAL PROPERTIES OF AXPS AND A POSSIBLE 
MODEL 

3.1 The diagram of spectral photon index and the accretion rate (F — M) 

In accretion-based models, a neutron star can be observed as an X-ray pulsar during 
a tracking phase. AXPs are supposed to rotate at a quasi-equilibrium period. And a 
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spin-down torque J is needed during this period. CHN suggests J changes with time 
according to the foUowing formula (Menou et al. 1999) 

j = m^ 2mRl^K{Rm) [l - „ "/p^ ] , (3) 

where VLK{Rm) is the Keplerian rotation rate at radius i?„j, and Rm is the magnetospheric 
radius which is determined by 

Rm ~ 6.6 X lO^sff m~2/7 

In Eq.21 m — M/A^Edd, where Msdd ~ 9.46 x lO^^gs"^ is the Eddington accretion rate, 
Bi2 = B/1012 Gauss. 

Following the arguments of Cannizzo et al. (1990), CHN suggests that after a dynam- 
ical time T the fossil disk loses mass self-similarly, 

mo , < t < r ; 

(5j 

mo{^) " , t>T ; 

where T ^ 10^'^ s, is the local dynamical time, and mo is a constant, which is normalized 
to the total initial disk mass, Md — Mddt, by mo = [(a — l)Md\/{aMEddT). a > 1 
is a constant that depends on directly the disk opacity (Francischell & Wijers 2003). 
Assuming an arbitrary value for a, equation [Sland|21can be combined to yield an analytic 
formula for Q{t), in terms of incomplete gamma functions (CHN). For a — 7/6, the 
solution of equation |51is 

f r!fc(i?,„,0)m^/', 0<<<T ; 
m={ ...3/7 (6) 

Since f]fe(i?„,0) = GAUs/Rt, then for i > T 

where Ci is a constant, and Ci = 25.48. If fl{t) evolves as equation[7| the spin period of 
the star nearly equals the evolving equilibrium period; this is a tracking phase (CHN). It 
is worth noting that the spin periods of SGRs and AXPs implied by equationd depend 
on both magnetic fields and mass accretion accretion rates. 

A strong propeller torque is needed to determine the spin evolution of the neutron star 
+ fossil disk system (CHN). However, the way in which a propelling neutron star loses 
angular momentum is not well understood (Francischell & Wijers 2003). In the most 
general form, the exact manner in which angular momentum is transferred between an 
accretion disk and a neutron star is a complex magneto-hydrodynamical problem, so there 
is no simple analytic solution for f2 changing with rh and B has been determined from 
equation 7. It is thought that most of the captured matter is ejected prior to reaching the 
neutron stars surface, i.e., in the strong propeller phase. However, after a sharp propeller 
cycle, the system reaches a tracking phase which is argued as producing "AXPs" phase 



anomalous X-ray pulsars 5 

(CHN). Clearly, employing different propeller torque (Francischell & Wijers 2003; Fabian 
1975; and CHN) will produce different results in determining the overall evolution of 
the system. How the propeller torques affect the period changing with time has been 
illustrated by Francischell & Wijers (2003). Based on above discussions, it is safe to 
assume that some kind of propeller torque acts on the neutron star between the strong 
propeller phase and tracking phase in order to spin it down with time. So, for simplicity, 
in this paper, we assume that the propeller torque is given by, 

j = m = -2imRl,Q.{t) . (8) 

where ^ is a constant, which is in the range o/ < ^ < 1 . ^ = 1 is for the strong propeller 
limit (Francischell & Wijers 2003), and ^ ^ is for the quasi- equilibrium phase, i.e., 
tracking phase. Then, from equatiorJHl we obtain 

n « CaBia^^™^^^ > (9) 

where C2 — 2.164 x 10^^*^^ is a constant. 

Combining equation[7|and|Hl and taking into account of il = 2tt/P and Cl = 2'kP / , 
the mass accretion rate as function of P and P is yielded 

7 

logm(F,F) = --logP + log P+ [10.26 -log ^] . (10) 
o 

Comparing equations 1101 with El a inverse relation between the photon index F and the 
mass accretion rate rn is inferred. Fitting the mass accretion rates m of 5 AXPs and 2 
SGRs estimated according to Eq.^|in the cases of ^ = 1, 0.5, and 0.05, respectively, we 
obtain 

logTO = ^-0.99r . (11) 

where A is a fitting constant, and A = 0.23, 0.53 and 1.53 corresponding to — \, 0.5 and 
0.05, respectively. Equation 11 shows a good correlation between F and m. The diagram 
ofV-m in the case of — 0.5 is plotted in Fig.|3| 

3.2 A Possible mechanism for the emission of AXPs 

Figure |3|implies that larger accretion rate produces the harder spectra (high-energy tail) 
in accretion-based model. The origin of this high-energy tail is unexplained very well at 
present. However, this high-energy tail could be due to thermal Comptonization by a hot 
coronal plasma, or it could be due to non-thermal emission (Tavani & Liang 1996). Based 
on the energetics of particle acceleration and cooling near the Alfven radius, Bottcher & 
Liang (2001) simulated the thermal-nonthermal radiation from a neutron star with both 
a weakly magnetized (~ lO^G) magnetospheric accretion shell and a normal magnetic 
case (~ lO^^G). The basic conclusion of Bottcher & Liang (2001) is: there is a positive 
correlation between spectral index and accretion rate, i.e., the spectral hardness and 
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accretion rate are anti-correlated. It is stressed that this conclusion is obtained in the 
situation of a neutron star with a weakly magnetic field, i.e., 10^ G. The interesting 
calculation for the photon spectra with the mass accretion rate is for the normal magnetic 
case (~ lO^^G), which are plotted in the figure 3 and figure 5 of Bottcher & Liang (2001). 
These two figures show that harder photon spectra are produced for higher accretion 
rates, particularly around energy lOkeV (2-10 keV, ACSA energy bands) in which we are 
concerned in this paper. As a natural consequence, in the case of a neutron star with 
a normal magnetic field, the photon index decreases with increasing accretion rates. In 
other word, there is a negative correlation between spectral index and accretion rates 
as observed in AXPs/SGRs. Furthermore, the hard X-ray spectral index resulted from 
the simulation of Bottcher & Liang (2001) are in excellent agreement with the values of 
F ~ 3 — 4 generally observed in AXPs. From this context, the fitting of F — m shows that 
the accretion rate is in the range from lO^'^ to 10^^, which is consistent with the model 
conditions suggested by Bottcher & Liang (2001). Therefore, we can apply the emission 
mechanism of Bottcher & Liang (2001) for the normal magnetic case to construct the 
accretion-powered emission model for AXPs, which can explain the relation between the 
accretion rate and spectral properties. 



4 DISCUSSIONS AND CONCLUSIONS 

The correlation between the spectra and the spin-down rate of AXPs/SGRs shows that 
the photon index decreases with increasing spin-down rate for AXPs and SGRs (see 
Fig. This fact has not so far been given a good explanation. If the increasing power- 
law emission with spin-down rate is consistent with the magnetar model, the spectra of 
AXPs/SGRs is expected to extend into the far-UV band (Marsden & White( 2001). This 
implies that observation of spectral breaks in the non-thermal persistent emission in the 
far-UV would be important evidence in support of the magnetar model. With respect 
to the accretion model (CHN), the spectral photon index and the accretion rate m (see 
Fig. 121 is investigated. To simplify the calculating of our model, an alternative propeller 
torque (see, Eq.(8)) to limit the spin-down rate at the boundary when the system transfers 
from the propeller phase to "tracking" phase is introduced. The analysis demonstrates that 
the correct spectral shape for the values of B and m can produce the rapid spin-down of 
AXPs and SGRs. A relation o/logm = 0.53 — 0.99r is derived for the case of ^ = 0.5. 
This relation plotted in figure 3 shows that the hardness of the spectra increases with 
increasing mass accretion rate. The same conclusion remains for different values of ^, 
only the slope of fitting lines is changed because ^ modifies the values of m (see Eg. (11)). 
In principle, the result shown in this paper is consistent with the model prediction given 
by Bottcher & Liang (2001) in the normal magnetic situation. Consequently, it is possible 
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that even stronger magnetic-field accreting neutron stars with low-mass accretion rate 
may be consistent with the data of AXPs/SGRs. 

Unlike the case for normal accretion onto neutron stars , the observed relationship 
between the absolute X-ray luminosity and the photon index F in AXPs/SGRs cannot 
be inferred from the relation of m — T. This is because the X-ray luminosity of the 
system is determined by the mass accretion rate on to the surface of the star, which 
could be different from the mass accretion rate through the disk (CHN). Two reasons 
can be involved to account for this. Firstly, in the propeller phase, it is not clear at all 
how much material eventually reaches the neutron star surface, and what is the radiation 
efficiency for the matter stopped the neutron star surface. Secondly, the X-ray luminosity 
may not cover the total luminosity. 

Acknowledgements We are particularly grateful to the anonymous referee for insightful com- 
ments and suggestions. We would like to thank Ling J.R. and Zhang W.M. of Tsinghua University 
for useful discussions for useful discussions. This work was supported by the National Natural 
Science Foundation of China under grant no. 10273011, the National 973 Project (NKBRSF 
G19990754), and the Special Funds for Major State Basic Research Projects. 



References 

Alpar, M.A., 2001, ApJ, 554, 1245 

Bottcher, M., & Liang, E.P., 2001, ApJ, 552, 248 

Cannizzo, J.K., Lee, H.M., & Goodman, J., 1990, ApJ, 351, 38 

Chatterjee, P., Hernquist, L., & Narayan, R., 2000, ApJ, 534, 373 (CHN) 

Corbet, R.H.D., & Mihara, T., 1997, ApJ, 475, L127 

Corbet, R.H.D.,et al., ApJ, 443, 786 

Dar, A., DeRujula, A., 2000, Results and Perspectives in Particle Physics, M., Greco ed.. Vol., 
XVII, 13 

Fabian, A.C., 1975, MNRAS, 173, 161 

FrancischeUi, G.J., & Wijers, R.A.M.J., 2003, astro-ph/0205212 

Gotthelf, E.V., & Vasisht, G., 1997, ApJ, 486, L133 

Gotthelf, E.V., Vasisht, G., & Dotani, T., 1999, ApJ, 522, L49 
Harding, A.K., et al. 1999, ApJ, 525, L125 
Heyl, J.S., & Hernquist, L., 1997, ApJ, 489, L67 
Hurley, K. et al., 1999, ApJ, 510, Llll 
Israel, G.L., et al., 1999, ApJ, 518, L107 

Kaspi, V.M., Chakrabarty, D., & Steinberger, J., 1999, ApJ, 525, L33 
Ho, W.C.G., & Lai, D., 2001, MNRAS, 304, L37 

Hurley, K., 2000, Proceedings 5th Huntsville GRB Symposium. AIP Conf. Series, 526, 763 

Lu, Y., & Cheng, K.S., 2002, ChJAA, 2, 161 

Marsden, D., et al. 2001, ApJ, 550, 397 

Marsden, D., & White, N.E., 2001, ApJ, 551, L155 

Mereghetti, S., & Stella, L., 1995, ApJ, 442, L17 

Menou, K., Esin, A.A., Narayan, R., Garcia, M.R., Lasota, J.-P., & McClintock, J.E., 1999, 
ApJ, 520, 276Me02 

Mereghetti, S., Chiarlone, L., Israel, G.L., & Stella, L. , 2002, |astro-ph/020512'2] 

Murakami, T. et al., 1999, ApJ, 510, L119 

Ozel, F., 2001, ApJ, 563, 276 

Paul, B., et al., 2000, ApJ, 537, 319 

Sonobe, T. et al., 1994, ApJ, 436, L23 

Sugizaki, M., et al. 1997, PASJ, 49, L25 

Tavani, M., & Liang, E.P., 1996, A&AS, 120, 133 

Thompson, C, & Duncan, R., 1996, ApJ, 473, 322 

Van Paradijs, J., Taam, R.E., & Van den Heuvel, E.P.J. , 1995, A&A, 299, L41 
White, N.E., Angelini, L., Ebisawa, K., Tanaka, Y., & Ghosh, P., 1996, ApJ, 463, L83 
Woods, P.M., et al., 1999, ApJ, 524, L55 
Woods, P.M., et al., 2000, ApJ, 535, L55 

Zane, S., TuroUa, R., Sella, L., & Treves, A., 2001, ApJ, 560, 384 



8 Y. Lu 

Table 1 SGR AND AXP TIMING PARAMETERS 



Object 


Start date 


p{sr 




r{BB + PLy 


References'* 


SGR 1900+14 


1998 Apr 30 


5.158971(7) 


61.0± 1.5 


2.1 


1, 2 




1998 Sep 16 


5.16025(2) 


61.0± 1.5 




2, 3 


SGR 1806-20 


1993 Oct 10 


7.468514(3) 


115.7± 0.2 


2.2 


4, 5 




1995 Oct 16 


7.46445(3) 


115.7+0.2 




4,5 


AXP 1048-59 


1994 Mar 03 


6.446646(1) 


32.9+0.3 


2.5 


6, 7 




1998 Jul 26 


6.45082(1) 


16.7 + 0.2 




7 


AXP 1841-05 


1993 Oct 11 


11.76668(6) 


41.3+0.1 


3.4 


8,9 




1998, Mar 27 


11.77243(7) 


41.3+0.1 




9 


AXP 2259+59 


1993 May 30 


6.97884(2) 


0.4883+0.0003 


3.9 


10,11 




1995, Aug 11 


6.9788793(8) 


0.4883+0.0003 




10 


AXP 0142+62 


1994 Sep 18 


8.68794(7) 


2.2+0.2 


4.0 


12 




1998 Aug 21 


8.68828(4) 


2.2+ 0.2 




7 


AXP 1709-40 


1996 Sep 03 


10.99758(6) 


22 + 6 


2.9 


13,14 



"Measured period (1 a error in last digit). 

''Assumed period derivative (from references). 

°BB+PL stands for black-body plus power-law spectral fit 

'* References.-(l)Hurely et al. 1999; (2)Woods et al. 1999; (3)Murakami et al. 1999; (4)Sonobe 
et al. 1994; (5)Woods et al. 2000; (6)Corbet & Mihara 1997; (7)Paul et al. 2000 ; (8)Gotthelf 
& Vasisht 1997; (9)Gotthelf, Vasisht, & Dotani 1999; (lO)Kaspi, Chakrabarty, & Steinberger 
1999; (ll)Cobet et al. 1995; (12)White et al. 1996; (13)Sugizaki et al. 1997; (14)Israel et al. 1999 
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Fig. 1 The T-Ct diagram of AXPs and SGRs. The triangles denote the AXPs 
and the circles are SGRs, and the symbols are all the same in the following 
figures. 
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Fig. 2 The r - B diagram of AXPs and SGRs. 




Fig. 3 The T — M diagram of AXPs and SGRs in the accretion model. The 
Une denotes the best fiiting relation of all the data points. 



